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The steady-state of heterogeneous catalysis, 
studied by first-principles statistical mechanics 
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The turn-over frequency of the catalytic oxidation of CO at RuO2(110) was calculated as function 
of temperature and partial pressures using ah mitio statistical mechanics. The underlying energetics 
of the gas-phase molecules, dissociation, adsorption, surface diffusion, surface chemical reactions, 
and desorption were obtained by all-electron density-functional theory. The resulting CO2 formation 
rate [in the full (T, pco,po2)-space], the movies displaying the atomic motion and reactions over 
times scales from picoseconds to seconds, and the statistical analyses provide insights into the 
concerted actions ruling heterogeneous catalysis and open thermodynamic systems in general. 

PACS numbers: 82.65.-(-r, 68.43.Bc, 68.43.De, 82.20.Db 



Under realistic conditions materials surfaces are in con- 
tact with a rich environment |Jj . Often, the resulting sur- 
face composition and surface-actuated material function 
are determined by equilibrium thermodynamics. A first- 
principles description is then possible with the "a& ini- 
tio atomistic thermodynamics" approach, that has been 
successfully applied to various systems (see e.g. Refs. 

IE IE S IB| and references therein). However, under 
many conditions equilibrium thermodynamics does not 
provide the appropriate description, and heterogeneous 
catalysis is a particularly interesting and important ex- 
ample. Here one is dealing with an open system, i.e., a 
supply of gases or liquids comes into contact with a solid 
surface, where a chemical reaction produces a new sub- 
stance that is then transported away. The entire concert 
of the various underlying, interlinked atomistic processes 
is in this case determined by kinetics. Still, the temper- 
ature and partial pressures of the reactants must be set 
such that the system runs under steady- state conditions. 
Only then the catalyst is not getting destroyed, but is sta- 
bly enhancing the rate of the desired chemical reaction. 
An ah initio description of the full steady-state situation 
of heterogeneous catalysis has not been achieved so far, 
and a microscopic understanding of the competing and 
concerting actions of the various atomistic processes is 
lacking. 

The present paper describes the "a6 initio statisti- 
cal mechanics" methodology appropriate for an open 
thermodynamic system with a continuous conversion of 
chemicals A and B into C, using the oxidation of CO 
at a RUO2 model catalyst as example. We employ 
density-functional theory (DFT) to obtain the energet- 
ics of all relevant processes, as there are: motion of the 
gas-phase molecules, dissociation, adsorption, surface dif- 
fusion, surface chemical reactions, and desorption. These 
calculations use the all electron, full-potential linear aug- 
mented plane wave (FP-LAPW) approach @, Q. The 
only notable approximations are the generalized gradi- 



ent approximation (GGA) for the exchange-correlation 
functional and the assumption that transition-state 
theory (TST) |^ is apphcable. In fact, both approxima- 
tions are well justified for the present study, and we will 
particularly address the GGA below when analyzing the 
results. A combination of DFT with TST, and subse- 
quently solving the statistical mechanicsproblem by the 
kinetic Monte Carlo (kMC) approach |sl ll(M. has been 
employed before (see, e.g. Refs. HOlIiEi). In 
distinction to (by now standard) "empirical kMC" cal- 
culations, that use just a few effective parameters, which 
have only limited (if any) microscopic meaning, these "a6 
initio kMC" calculations include an extended set of el- 
ementary processes with full and direct physical mean- 
ing. And describing an open system with a continuous 
conversion of chemicals A and B into C implies addi- 
tional complexity. For the example discussed below (the 
catalytic oxidation of CO) we will solve the statistical 
mechanics problem, gaining microscopic insight into the 
full dynamics from picoseconds to microseconds and even 
to seconds. The reported results demonstrate the new 
quality of and the novel insights gained by such descrip- 
tion. We will also compare the results to those obtained 
by "constrained thermodynamics" 3]. Although (as ex- 
pected) noticeable differences occur for certain environ- 
mental conditions, we will confirm that the much simpler 
"constrained thermodynamics" approach can provide im- 
portant guidance on where in (T, p)-space high reaction 
rates are likely to be expected, and when results obtained 
for certain (T,p) conditions can be extrapolated to oth- 
ers. 

The now obtained turn-over-frequencies (TOFs) of the 
catalytic oxidation of CO are in unexpected and unprece- 
dented agreement with experimental results of Peden and 
Goodman Ti^ and Wang et al. . This is found to be a 
consequence of the fact that the high reaction conditions 
are not just ruled by a singular chemical reaction path- 
way. Instead, for the steady-state and high TOF condi- 
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tions it is necessary that various processes play together 
in a most efficient manner. For the present example one 
crucial point is to realize an optimum disordered and dy- 
namic "phase" at the surface. Modest errors due to the 
approximate treatment of the exchange-correlation func- 
tional (e.g. the GGA) are then not crucial, as long as 
the trends of the energetics of the various interplaying 
processes are described correctly. In this sense, the fre- 
quently requested chemical accuracy for the description 
of individual processes appears to be a misleading con- 
cept. At least for the present system, a careful combina- 
tion of DFT and statistical mechanics is more important. 

The methodology will be used to study CO oxidation 
over RuO2(110) (the surface is sketched in Fig. QJi), as 
this system has recently received considerable attention 
as a highly active model catalyst (see Ref. and refer- 
ences therein). In fact, this was previously called the Ru 
catalyst, but recent experimental and theoretical work 
has shown that at realistic O2 pressure the Ru(OOOl) sur- 
face is transformed into an epitaxial RuO2(110) film (see 
Ref. and references therein). It is by now also estab- 
lished that this RuO2(110) surface actuates the catalytic 
reaction, and, although domain boundaries and steps are 
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FIG. 1: a) Perspective view of the RuO2(110) surface, illus- 
trating the two prominent adsorption sites in the rectangular 
surface unit cell. These sites are labeled as br (bridge) and cus 
(coordinatively unsaturated) site, and both are occupied in 
the exarnple with oxygen atoms. This is the "high po^ termi- 
nation" b) Time evolution of the site-occupation num- 
bers from the thermodynamic surface termination shown in 
a) to the steady-state of catalysis [cf. snapshot c)] . The tem- 
perature is r = 600 K and the CO and O2 partial pressures 
are that of the optimum reaction conditions: pco ~ 20 atm, 
PO2 = 1 atm. c) Snapshot (top view of the RuO2(110) sur- 
face) from movies displaying the dynamics of the surface un- 
der realistic catalytic conditions. Though the full atomic 
structure is considered in the calculations, for clarity we have 
marked here the substrate bridge sites only by gray stripes 
and the cus sites by white stripes. Oxygen adatoms are drawn 
as red circles and adsorbed CO molecules as blue circles. 
Movies for various pressure conditions for runtimes of 500 ns, 
as well as 1 s can be found at http ://w3. rz-berlin.mpg.de/ 
~reut er /movies / movies . ht ml . 



TABLE I: DFT binding energies, Eb, for CO and O (with 
respect to (l/2)02) at br and cus sites (cf. Fig. ^), diffusion 
energy barriers, AE^^g, to neighboring br and cus sites, and 
reaction energy barriers, AiSj'oac) of neighboring species at br 
and cus sites. All values are in eV. 
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present, their influence is not significant [17|. The surface 
unit cell is rectangular and contains two adsorption sites 
[sl : bridge (br) and coordinatively unsaturated (cus) sites 
(cf. Fig. ^). Either of them can be empty or occupied 
by O or CO, and adsorbate diffusion can go br-to-br, br- 
to-cus, cus-to-cus, or cus-to-br. Since this comprises the 
possibility of missing O*"^ atoms, we note that our treat- 
ment implicitly includes the effect of O surface vacancies. 
A total of 26 different elementary processes are possible 
on this lattice, and all were carefully analyzed by DFT to 
obtain their pathways and energy barriers 0;lla|- TableU 
summarizes the adsorption energies and the diffusion and 
reaction energy barriers used for the kMC study, as ob- 
tained from DFT calculations with a (1x2) surface unit- 
cell. In a systematic study of CO and O (co)adsorption 
at various coverages we found adsorbate-adsorbate inter- 
action to be always smaller than 150 meV. Thus, lateral 
interactions in this system are small (compared to the 
other energies), and will therefore be neglected. 

CO adsorption into vacant cus or bridge sites is non- 
dissociative, while oxygen adsorption is dissociative and 
requires two vacant neighboring sites, i.e. a br-br, cus- 
cus, or br-cus pair. The adsorption rate per free site 
is given by the local sticking coefficient, S, and the ki- 
netic impingement: Fad = S p/ . Here m is 
the mass of the gas-phase molecule, and k-Q is the Boltz- 
mann constant. The rate of the time reversed process 
(desorption) then follows from the relation: Fdes/Fad = 
exp {{Fh — fi)/{kBT)), where Fb is the free energy of the 
adsorbed species (approximated by Eb), and ^l^^p) is 
the chemical potential of the gas-phase molecule . The 
S(T) are thus obtained from the calculated total en- 
ergy surfaces of desorption together with detailed bal- 
ance. The only uncertainty here arises from the vibra- 
tional properties of the transition state which translates 
into uncertainties in the desorption rate by a factor of 10 
(at most 100). For the diffusion processes we use a prcf- 
actor of 10"'^ Hz, which has an uncertainty of a factor of 
10. We carefully checked that these uncertainties do not 
affect our below reported results and conclusions. Details 
of the employed new methodology and of the various test 
calculations will be published elsewhere Detailed 
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balance of the scenario was carefully checked by confirm- 
ing that the earlier results obtained by "atomistic ther- 
modynamics" Q are exactly reproduced by the present 
statistical mechanics treatment, if surface reactions are 
not allowed to occur. 

Kinetic Monte Carlo runs were performed for about 
1000 different (T, PC0)P02) conditions covering the tem- 
perature range 300 K < T < 800 K and partial pressures 
from 10~^° to 10'^ atmospheres. Several calculations were 
done on a system with (30 x 30) surface sites (450 bridge 
plus 450 cus sites), but the vast amount of calculations 
was done for a 400 sites system. The results for both 
system sizes were identical. The kMC simulations were 
run until steady-state is reached (cf. Fig. Then 
the movies were recorded (cf. Fig. ^) and the statistical 
analyses of the frequencies of the various elementary pro- 
cesses performed. The latter also gives the total TOFs 
of the CO2 formation. 

The results show that catalytic reaction conditions are 
only established when both CO and O2 partial pressures 
exceed certain values below which the surface is in a 
thermodynamic equilibrium phase of RuO2(110), charac- 
teristic for low oxygen pressures and routinely observed 
under ultra-high vacuum (UHV) conditions: All bridge 
sites are occupied by oxygen atoms and all cus sites are 
empty. For higher pressures three surface conditions are 
worth mentioning, which are e.g. obtained for T = 600 K, 
PO2 — 1 atm, and pco — 10~^, 20, and 10"^ atm. 

(i) At the low pco, all bridge and all cus sites are 
covered with oxygen atoms (cf. Fig. QJi). This is essen- 
tially the thermodynamic high po^ phase 0, 0| ■ There is 
a noticeable desorption/adsorption dynamics, i.e., about 
every 40 /zs two of the 900 O adatoms of the (30 x 30) 
simulation cell desorb (mainly from cus sites) as O2. The 
resulting vacancies are then rapidly filled again (within 
1 ns), most of the time with oxygen. Only rarely CO ad- 
sorbs, and even if it does, it rather desorbs again than 
initiating a reaction. The overall CO2 formation rate un- 
der these conditions is with 0.9-10"'^^ cm^^s^^ very low. 

(a) For Pco = 20 atm the situation corresponds to 
what was previously suggested to be that of high chemi- 
cal activity , and this suggestion is now confirmed. The 
time to reach steady-state is remarkably long, in partic- 
ular when compared to the pico second timescale of the 
underlying atomistic processes (cf. Fig. QJd). Finally, an 
interesting mix of CO molecules and O atoms at bridge 
and cus sites is established (a typical snapshot is shown in 
Fig- Gt), with average occupation numbers Nqq = 0.11, 
iVg'g = 0.70, iV^'' = 0.89, and A^g"'' = 0.29. The dy- 
namics of this surface is extremely fast. It is mainly due 
to CO desorption and adsorption. The rate of CO2 for- 
mation is a factor of 0.0004 lower, but truly significant, 
namely 4.6 • lO^^cm"^ s^^. It is dominated by the reac- 
tion 00'=^'' + O™" — > CO2. The CO^^' + C^"" — > CO2 
and CO'^"^ -|- O'"' — > CO2 reactions also contribute, 
though their rates are by a factor of 0.30 and 0.01 lower 



than the first one. The CO'"' -I- O^'' — > CO2 reaction 
is insignificant. This result is different from what one 
would expect from the reaction energy barriers (cf. table 
m, which would give the lead to CO'''' -I- O'^"'', demon- 
strating the importance of the proper mix of surface-site 
occupations. 

(in) When the CO pressure is increased further, to 
Pco = 10^ atm, the surface becomes fully covered with 
CO. This result is at variance with the "constrained ther- 
modynamics" study in which the CO2 formation was for- 
bidden. Now we find that catalysis is practically poisoned 
by adsorbed CO. The encountered situation is in fact 
close to that where the RUO2 catalyst will be reduced by 
CO to the pure Ru metal. 

Figure 13 summarizes the results by showing the var- 
ious steady-state surface structures as well as a map of 
the TOFs in {pco, po2)-space. The highest activity is 
found to be in very good agreement with the early ex- 
perimental results by Peden and Goodman and oc- 
curs whenever the environmental conditions lead to the 
dynamic coexistence " pha se" described above under (ii). 
Recently, Wang et al. [l5j have also measured the TOFs 
for the RuO2(110) surface at T = 350K for various pres- 
sures, and in Fig. |31we compare their results with ours. 
The agreement is again far better than what one would 
have expected: The theoretical and experimental TOF 
values at the optimum pressures are practically identi- 
cal, and also the optimum pressure conditions (when the 
"dynamic phase" is realized) agree very well. Because the 
errors due to the GGA are in the range of 0.1-0.3 eV, this 
good agreement between theory and experiment, seen in 
Fig. 13 may seem fortuitous. However, it is worth noting 
that the position of the optimum catalytic efficiency in 
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FIG. 2: a) Steady-state surface structures obtained by ab 
initio kMC calculations at T = 600 K and various CO and 
O2 partial pressures. In all non- white areas, the average site 
occupation is dominated (> 90%) by one species, i.e. ei- 
ther O, CO or empty sites (— ). b) Map of the corresponding 
turn-over frequencies (TOFs) in cm~'^s~^: White areas have 
a TOF < 10^^cm~^s~^, and each increasing gray level repre- 
sents one order of magnitude higher activity. Thus, the dark- 
est region corresponds to a TOF higher than 10^^cm~^s~^. 
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FIG. 3: Rate of CO2 formation at T = 350 K. The experi- 
mental steady-state results of Wang et al. [irI | are presented 
as dotted lines, and the theoretical results are shown as green, 
solid lines. Rates given as function of PO2 atpco = 10"^" atm 
(left), and as function of pco at P02 = 10~^° atm (right). 



(T, p)-space and to some extent also the value of the TOF 
are not determined by the energetics of a singular process 
alone, but by the action of many players. Apparently, the 
DFT-GGA calculations describe the differences between 
the various surface processes better than the individual 
absolute values. In this respect it is important to realize 
that a combination of different calculations (employing 
different approximations) or of theory and experiment 
could have spoiled the description. The consistent treat- 
ment of all participating processes implies that the opti- 
mum mix of O and CO at the surface is described well: 
The abundance of C0'^"*'-0'^"'' nearest neighbor pairs (as 
well as C0'"'-0^"** nearest neighbor pairs) is apparently 
playing a role of similar importance as the energy barri- 
ers. Particularly at the optimum TOF conditions there 
is also an effective compensation, e.g. when too high ad- 
sorption energies result in enhanced adsorption, but also 
in reduced reaction barriers etc. Away from the opti- 
mum TOF conditions such compensation effects become 
less effective, and modest differences between the theo- 
retical and experimental results arise. Here DFT-GGA 
errors are more influential, and for the experimental data 
we expect that contributions from surface imperfections 
may play a bigger role. 

In summary, we computed the surface kinetics of CO 
oxidation catalysis at RuO2(110). The TOFs are pre- 
sented in (r,pcOjP02)"Space, clearly identifying a nar- 
row region of highest catalytic activity. In this region ki- 
netics builds an adsorbate composition that is not found 
anywhere in the thermodynamic surface phase diagram. 
The statistical analysis of the surface dynamics and of the 
various processes reveals several surprising results. For 
example, the chemical reaction with the most favorable 
energy barrier happens a factor of 0.30 less frequently 
than the energetically second favorable reaction. 

The results also clarify how and when a bridging of the 
pressure gap between UHV studies and realistic pressure 
conditions is possible. The considered system has in fact 
two important components to this issue: At first, an O2- 



rich environment changes the material from Ru to RUO2. 
Thus, earlier high pressure studies on the Ru catalyst 
were actually looking at RUO2. Second, after Ru02 has 
been formed it is important to know how to set the pres- 
sure conditions correctly, as also on RuO2(110) there is a 
low-pressure surface phase (the orange region in Fig. |2t), 
that has little in common with the catalytically active sit- 
uation. The obtained agreement between the theoretical 
results and experimental data confirms furthermore (cf. 
also Ref. that CO2 is primarily formed from ad- 

sorbed CO and O, and that the metal-oxide, once it was 
created, does not play an active role, i.e., there is no in- 
dication of significant bulk diffusion. Thus, the catalysis 
is explained in terms of a Langmuir-Hinschlwood mecha- 
nism. Some contributions via the Eley-Rideal mechanism 
(e.g. scattering of gas-phase CO at O™^) may play a mi- 
nor role, but this process requires further calculations and 
experiments. With the advent of first-principles TOF- 
maps obtained from DFT calculations (cf. Fig. EJ, a 
more detailed data base is in general becoming available 
for comparison with experiments, which will eventually 
advance the microscopic understanding of catalysis. 
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